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Abstract. The first sources of ionizing radiation to condense out of the 
fs^ I dark and neutral IGM sent ionization fronts sweeping outward through 

their surroundings, overtaking other condensed objects and photoevapo- 
rating them. This feedback of universal reionization on cosmic structure 
formation is demonstrated here by gas dynamical simulations, including 
radiative transfer, for the case of a cosmological minihalo of dark matter 
Tij- . and baryons exposed to an external source of ionizing starlight, just after 

I the passage of the global ionization front created by the source. 

o 
o 

2 . The neutral, opaque IGM out of which the first bound objects condensed was 

I dramatically reheated and reionized at some time between a redshift z ~ 50 

Ci ■ and z ~ 5 by the radiation released by some of these objects. When the first 

sources turned on, they ionized their surroundings by propagating weak, R-type 
- ionization fronts which moved outward supersonically with respect to both the 

/\ I neutral gas ahead of and the ionized gas behind the front, racing ahead of the 

■ hydrodynamical response of the IGM, as first described by Shapiro (1986) and 

Shapiro & Giroux (1987). These authors solved the problem of the time- varying 
radius of a spherical I-front which surrounds isolated sources in a cosmologically- 
expanding IGM analytically, taking proper account of the I-front jump condition 
generalized to cosmological conditions. They applied these solutions to deter- 
mine when the I-fronts surrounding isolated sources would grow to overlap and, 
thereby, complete the reionization of the universe (Donahue & Shull 1987 and 
Meiksen &: Madau 1993 subsequently adopted a similar approach to answer that 
question). The effect of density inhomogeneity on the rate of I-front propagation 
was described by a mean "clumping factor" q > 1, which slowed the I-fronts by 
increasing the average recombination rate per H atom inside clumps. This suf- 
fices to describe the rate of I-front propagation as long as the clumps are either 
not self-shielding or, if so, only absorb a fraction of the ionizing photons emitted 
by the central source. Numerical radiative transfer methods are currently under 
development to solve this problem in 3D for the inhomogeneous density distribu- 
tion which arises as cosmic structure forms, so far limited to a fixed density field 
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Figure 1. (a) (left) MINIHALO INITIAL CONDITIONS BEFORE 
REIONIZATION: (Top) gas density; (Bottom) gas velocity versus dis- 
tance from minihalo center, (b) (right) I-FRONT PHOTOEVAPO- 
RATES MINIHALO: (Top) I- front position along x-axis versus time; 
(Bottom) Mass fraction of the initial mass Mi of minihalo hydrostatic 
core which remains neutral (H I) versus time. 

without gas dynamics (e.g. Abel, Norman, & Madau 1999; Razoumov & Scott 
1999). The question of what dynamical effect the I- front had on the density 
inhomogeneity it encountered, however, requires further analysis. 

2. The Photoevaporation of Dwarf Galaxy Minihalos Overtaken by 
a Cosmological Ionization Front 

Shapiro and Raga (2000) describe radiation-hydrodynamical simulations of the 
photoevaporation of a cosmological minihalo overrun by a weak, R-typc I-front 
in the surrounding IGM, created by a quasar with emission spectrum Fi, oc u^^'^ 
{v > i'h) and a flux at the location of the minihalo which, if unattenuated, would 
correspond to A^phjSe/'^Mpc ~ where rMpc is the distance (in Mpc) between 
source and minihalo and iVph,56 is the H-ionizing photon luminosity (in units of 
10^^ s~^). For further discussion and references to previous work, the reader is 
referred to this companion paper by Shapiro and Raga. Our simulations in 2D, 
axisymmetry used an Eulerian hydro code with Adaptive Mesh Refinement and 
the Van Leer flux-splitting algorithm, which solved nonequilibrium ionization 
rate equations (for H, He, C, N, O, Ne, and S) and included an explicit treatment 
of radiative transfer by taking into account the bound-free opacity of H and He. 
The heavy element abundance was 10~^ times solar. 

Here we present new results of identical simulations, except with a stellar 
source instead of a quasar-like source, with a 50,000 K blackbody spectrum and 
with luminosity and distance adjusted to keep the ionizing photon flux the same 
(i.e. A'ph^se/j^Mpc ~ -'-)• initial condition before ionization, shown in Figure 
1(a), is that of a lO^M© minihalo in an Einstein-de Sitter universe (Ocdm = 
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Figure 2. PHOTOEVAPORATING MINIHALO. 50 Myr after turn- 
on of stellar UV source to the left of computational box along the x- 
axis. (a) (Upper Panel) isocontours of atomic density, logarithmically 
spaced, in (r,x)— plane of cylindrical coordinates; (b) (Lower Panel) 
velocity arrows are plotted with length proportional to gas velocity. 
An arrow of length equal to the spacing between arrows has velocity 
25kms^^. Solid line shows current extent of gas originally in hydro- 
static core. Dashed line is I- front (50% H- ionization contour). 



1 — ^bary; ^bary^^ = 0.02; h = 0.7) which collapses out and virializes at ZcoU = 9, 
yielding a truncated, nonsingular isothermal sphere of radius Rc = 0.5 kpc in 
hydrostatic equilibrium with virial temperature Tvir = 5900 K and dark matter 
velocity dispersion a = 6.3kms~^, according to the solution of Shapiro, Iliev, & 
Raga (1999), for which the finite central density inside a radius about 1/30 of the 
total size of the sphere is 514 times the surface density. This hydrostatic core of 
radius Rc is embedded in a self-similar, spherical, cosmological infall according 
to Bertschinger (1985). 

The results of our simulation on an (r, x)-grid with 256 x 512 cells (fully 
refined) are summarized in Figures 1-5. The background IGM and infalling 
gas outside the minihalo [centered at {r,x) = (0,7.125 x 10^^ cm)] are quickly 
ionized, and the resulting pressure gradient in the infall region converts the infall 
into an outflow. As expected, the hydrostatic core of the minihalo shields itself 
against ionizing photons, trapping the I- front which enters the halo, causing it to 
decelerate inside the halo to the sound speed of the ionized gas before it can exit 
the other side, thereby transforming itself into a weak, D-type front preceded by 
a shock. The side facing the source expels a supersonic wind backwards towards 
the source, which shocks the IGM outside the minihalo, while the remaining 
neutral halo material is accelerated away from the source by the so-called "rocket 
effect" as the halo photoevap orates (cf. Spitzer 1978). Since this is a case of 
gas bound to a dark halo with a < lOkms"^, this photoevaporation proceeds 
unimpeded by gravity. Figure 1(b) shows the position of the I- front inside the 
minihalo as it slows from weak, R-type to weak D-type led by a shock as it 
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Figure 3. PHOTOEVAPORATING MINIHALO. One time-slice 50 
Myr after turn-on of stellar UV source located to the left of computa- 
tional box along the x-axis. Prom top to bottom: (a) isocontours of 
pressure, logarithmically spaced, in (r,x)— plane of cylindrical coordi- 
nates; (b) pressure along the r = symmetry axis; (c) temperature; (d) 
H I fraction; (e) He I (solid) and He II (dashed) fractions; (f) bound-free 
optical depth along r = axis at the threshold ionization energies for 
H I (solid) , He I (dashed) , He II (dotted) . Key features of the flow are 
indicated by the numbers which label them on the temperature plot: 
1 = IGM shock; 2 = contact discontinuity between shocked halo wind 
and swept-up IGM; 3 = wind shock; between 3 and 4 = supersonic 
wind; 41 = I-front; 4S = shock which precedes I-front; 5 = boundary of 
gas originally in hydrostatic core; 6 = shock in shadow region caused 
by compression of shadow gas by shock-heated gas outside shadow. 
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Figure 4. OBSERVATIONAL DIAGNOSTICS OF PHOTOEVAP- 
ORATING MINIHALO I: IONIZATION STRUCTURE OF METALS. 
Carbon, nitrogen, and oxygen ionic fractions along symmetry axis at 
t = 50 Myr. 



advances across the original hydrostatic core. Figure 1(b) also shows the mass 
of the neutral zone within the original hydrostatic core shrinking as the minihalo 
photoevaporates within about 100 Myrs. This photoevaporation time is 50 % 
larger than that found previously for the quasar case. Figures 2 and 3 show 
the structure of the photoevaporative flow 50 Myrs after the global I-front first 
overtakes the minihalo, with key features of the flow indicated by the labels on 
the temperature plot in Figure 3. A strong shock labelled "4S" clearly leads 
the D-type I-front (labelled "41") as it advances through the minihalo core, by 
contrast with the quasar case in which hard photons penetrated deeper into 
the neutral gas and preheated it, thereby weakening the shock which leads the 
I-front. The softer stellar spectrum also explains why helium on the ionized 
side of the I-front is mostly He II, rather than He III as in the quasar case, 
while the neutral side is completely He I, rather than a mix of He I and II as in 
that case. Figure 4 shows the spatial variation of the relative ionic abundances 
of C, N, O ions along the symmetry axis after 50 Myrs. While the quasar 
case showed the presence at 50 Myrs of low as well as high ionization stages 
for the metals, the softer spectrum of the stellar case yields less highly ionized 
gas on the ionized side of the I-front (e.g. mostly C III, N III, O III) and the 
neutral side as well (e.g. C II, N I. O I and II). The column densities of H I, 
He I and II, and C IV for minihalo gas of different velocities as seen along 
the symmetry axis at different times are shown in Figure 5. At early times, 
the minihalo gas resembles a weak Damped Lyman a ("DLA") absorber with 
small velocity width {'^ lOkms^^) and A^hi ^ lO^'^cm"^, with a Lyman-a- 
Forest("LF")-like red wing (velocity width ^ lOkms"^) with A/ri ^ lO^^cm"^ 
on the side moving toward the source, with a He I profile which mimics that 
of H I but with A^Hei/A^Hi ~ [He]/[H], and with a weak C IV feature with 
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Figure 5. OBSERVATIONAL DIAGNOSTICS OF PHOTOEVAP- 
ORATING MINIHALO II. ABSORPTION LINES: Minihalo column 
densities (cm~^) along symmetry axis at different velocities. (Top) 
H I; (Middle) He I (solid) and He II (dotted); (Bottom) C IV. Each 
box labelled with time (in Myrs) since source turn-on. 



New ~ lO^^cm"^ displaced in this same asymmetric way from the velocity of 
peak H I column density. After 160 Myr, however, only a narrow H I feature 
with LF-like column density A^hi ~ lO^^cm"^ remains, with Nuei/Nm ~ 1/4, 
A^Heii/A^Hi ~ 10^, and Ncw/Nni ~ 3[C]/[C]0. Future work will extend this 
study to minihalos of higher virial temperatures, for which gravity competes 
more effectively with photoevaporation. 
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